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The structure of shocks that form at the exhaust boundaries during collisionless reconnection 
of anti-parallel fields is studied using particle-in-cell (PIC) simulations and modeling based on the 
anisotropic magnetohydrodynamic equations. Large-scale PIC simulations of reconnection and com- 
panion Riemann simulations of shock development demonstrate that the pressure anisotropy pro- 
duced by counterstreaming ions within the exhaust prevents the development of classical Petschek 
switch-off-slow shocks (SSS). The shock structure that does develop is controlled by the firehose 
stability parameter e = 1 — /j.o(P\\ — P±)/B'^ through its influence on the speed order of the in- 
termediate and slow waves. Here P\\ and P± are the pressure parallel and perpendicular to the 
local magnetic field. The exhaust boundary is made up of a series of two shocks and a rotational 
wave. The first shock takes e from unity upstream to a plateau of 0.25 downstream. The condition 
e = 0.25 is special because at this value the speeds of nonlinear slow and intermediate waves are 
degenerate. The second slow shock leaves s = 0.25 unchanged but further reduces the amplitude of 
the reconnecting magnetic field. Finally, in the core of the exhaust e drops further and the transition 
is completed by a rotation of the reconnecting field into the out-of-plane direction. The acceleration 
of the exhaust takes place across the two slow shocks but not during the final rotation. The result is 
that the outflow speed falls below that expected from the Walen condition based on the asymptotic 
magnetic field. A simple analytic expression is given for the critical value of e within the exhaust 
below which SSSs no longer bound the reconnection outflow. 

I. INTRODUCTION 

Magnetic reconnection plays an important role in converting magnetic energy to plasma kinetic and thermal energy. 
The bulk of the energy released during reconnection takes place downstream of the magnetic X-line where reconnected 
field lines expand outwards to release their tension. Further, the geometry of the outflow exhaust ultimately determines 
the rate of reconnection, the narrow exhaust of the Sweet-Parker model producing slow reconnection and the open 
outflow exhaust proposed by Petschek producing fast reconnection [1 . Petschek specifically proposed that a pair of 
switch-off slow shocks (SSS) stand in the inflow into the exhaust and convert the released magnetic energy into the 
Alfvenic outflow and thermal energy. Magnetohydrodynamic simulations of reconnection with a localized resistivity 
imposed at the X-line confirmed the development of Petschek's open outflow configuration and the formation of SSSs 
at the exhaust boundaries \2^. Particle heating by these shocks has been proposed to cause the X-ray emission in 
solar flares |4] . 

Magnetic reconnection also takes place in environments in which collisions are either weak or essentially absent. 
Observations of reconnection in the solar wind in particular seem to suggest that reconnection X-lines and associated 
exhausts grow to very large scales and resemble the open outflow geometry predicted by Petschek [5]. On the other 
hand, direct observations of SSSs in the Earth's magnetosphere and the solar wind are infrequent [6]. 

An important result of the simulations of collisionless reconnection and the resultant Hall reconnection model is the 
spontaneous formation of an open outflow configuration analogous to that predicted by Petschek. This configuration 
develops even without an adhoc localized resistivity [3 H] . A key question is whether the fast rates of reconnection 
seen in the relatively small systems explored in simulations persist in very large systems. Scaling studies seem to 
suggest that collisionless reconnection remains fast in very large systems (9H11). On the other hand, while Hall MHD 
simulations of collisionless reconnection reveal the formation of SSSs at the boundaries of reconnection exhausts, 
kinetic simulations (hybrid or PIC models) reveal that the reconnecting magnetic field never switches off as expected 
from the SSS model [12 . Since it is the release of magnetic energy downstream from the X-line that ultimately 
drives the outflow rather than the dynamics close to the X-line, the absence of the SSS in kinetic simulations calls 
into question the conjecture that fast collisionless reconnection actually can scale to very large systems. Thus, a key 
requirement for demonstrating that Hall reconnection can explain the fast release of energy that takes place in large 
systems is to pin down the specific mechanism driving the Alfvenic outflow. 

It is well known from observations [TsHKj and modeling HSl Hi] of anti-parallel magnetic reconnection that ions 
counterstream through the exhaust with relative velocities of the order of the Alfven speed, increasing the plasma 
pressure parallel to the local magnetic field. The consequence of this increase in parallel pressure is the development 
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of a barrier in the pseudo-potential that describes the slow shock transition. This prevents the slow shock from 
switching-off the reconnecting magnetic field [18l [19] . For highly oblique slow shocks relevant to reconnection and 
where the upstream value of e is unity the critical value of the firehose stability parameter e^r below which the SSS 
solution becomes inaccessible is given by 

where is the ratio of the plasma to magnetic pressure upstream. If e falls below this critical value anywhere within 
the exhaust, the SSS solution becomes inaccessible, e does fall below this value in kinetic reconnection simulations 
[18]. The presence of the potential barrier to the formation of the SSS actually arises because of the formation of a new 
slow shock solution, which is characterized by intermediate Mach number M/ = ^JV^ii^pj sB^ = 1 both upstream 
and downstream [19 . We denote this as a 1,1 slow shock or 11-SS. This 11-SS can also switch off the magnetic 
field downstream but differs from the SSS where the downstream intermediate Mach number Mjd < 1. However, we 
show in the present manuscript that this shock solution cannot reduce e below the critical value given in Eq. ([T]) and 
therefore can not explain the values of s below Scr seen in the reconnection simulations p!8|. 

Riemann simulations of the structure of slow shocks which a constant normal magnetic field is added to 

a Harris equilibrium [20l |2T] , have revealed that the value of £ downstream of the shock transition tends to form 
a plateau at a value of 0.25 and it was demonstrated that at this value the phase speeds of nonlinear intermediate 
and slow waves are degenerate [19 . It was suggested that the coupling to a rotational mode at e = 0.25 produced 
turbulence at the ion inertial scale di = c/ujpi that produced sufficient scattering of ions to prevent e from falling 
below 0.25. However, neither the SSS nor the 11-SS can reduce e from unity upstream to 0.25 downstream. 

Here we present the results of a PIC reconnection simulation that is large enough to reveal the formation of an 
e = 0.25 plateau. The structure of the shocks that make up reconnection exhausts are explored by comparing 
the results of these reconnection simulations with parallel Riemann simulations and analytic analysis based on the 
anisotropic MHD equations. We show that the exhaust boundary is defined by two slow shock transitions followed 
by an RD. An anomalous slow shock (A-SS) [22], in which Mj > 1 both upstream and downstream, reduces the 
upstream tangential magnetic field and takes e from unity upstream to 0.25. A second slow shock forms on the 
s = 0.25 plateau and further reduces the tangential magnetic field. This degenerate slow shock (D-SS) is a special 
case of the 11-SS corresponding to Mj = Msl = 1 both upstream and downstream, where Msl is the slow mode 
Mach number. The D-SS forms a compound wave with a RD whose rotation onsets when e begins to fall below 0.25. 
In the reconnection simulation the core of the exhaust becomes firehose unstable, which prevents further acceleration 
of the exhaust outfiow by the RD. As a consequence, the exhaust velocity falls below the value expected based on the 
upstream Alfven speed. 



II. SIMULATION DETAILS AND RESULTS 



The PIC simulation is performed with the PSD code [23] with periodic boundary conditions. The initial state is a 
double-Harris configuration. However we focus on one current sheet of = —Bat8inh{x/wi) and n = nhsech^{x/wi)-\- 
Ua- The magnetic field is normalized to the asymptotic magnetic field Ba^ the density to the Harris density n/^, 
velocities to the Alfven speed Ca = Ba/ y^/io^i^^, lengths to the ion inertial length di = ^J rrii / (iQifihe^ ^ times to 
the inverse ion cyclotron frequency rt~^ = rrii/Bae^ and temperatures to miC\. Other important parameters are 
mi/rrie = 25, c = 15, = 0.2, Wi = 2 and uniform initial T^^g = 0.25, which imply that Pa = 0-2. The system size is 
S19.2di X 409. 2(i^ resolved by grids 16384 x 8192 with 100 particles per cell. Particles are advanced with a time step 
At = 0.01. 

In Fig. [TJa) is the out-of-plane current density Jy from a simulation with quasi-steady reconnection. The X-line 
is at the origin. The black contours are the in-plane magnetic field lines. The outfiow exhaust forms a Petschek-like 
open-outfiow configuration. In (b) we show the firehose stability parameter s. The counterstreaming ions that have 
been widely documented in satellite observations of reconnection exhausts p^HT5] drive the anisotropy toward the 
firehose threshold. The blue curves bound the region where e < 1. The red color table marks the region where s falls 
between 0.1 and 0.4. In the central region of the exhaust, marked by the green contour, is the firehose unstable region 
where e < 0.0. The significant structure in the firehose unstable region suggests that this instability is active in the 
core of the exhaust |T2l [24]. With increasing distance from the X-line the region where 0.1 < £ < 0.4 widens. In 
(c) we show cuts of s across the exhaust at four locations as marked in (b). The tendency to form a plateau at 0.25 
is evident in spite of the turbulence from the firehose unstable region. Petschek's reconnection model requires the 
development of a standing SSS that propagates in the infiow frame at the phase speed of intermediate wave. Since 
the intermediate wave speed goes to zero when e goes to zero, the formation of a slow shock at the exhaust boundary 
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is going to differ substantially from the traditional SSS. Similar plateaus at £ = 0.25 have been identified in Riemann 
simulations of slow shocks [18 . 

The structure of the exhaust and its boundaries are shown in greater detail in Fig. [2] where various physical 
quantities are plotted along the black cut in Fig. [lib). The gradual drop in Bz (in (b)) beginning at the exhaust 
boundary and the corresponding rise in pressure (in (d)), with the total plasma and magnetic pressure (in (d)) 
being nearly constant, is clear evidence for the development of a slow shock at the exhaust boundary. However, the 
downstream does not appear to be switched-off (i.e., B^ does not drop to zero and remain zero in the core of the 
exhaust) as in a MHD SSS. The out-of-plane component By indicates the presence of the residual of the right-handed 
polarized dispersive whistler wave that coincides with the front of the SS transition. These whistler waves open the 
reconnection dissipation region into an Petschek-like configuration [25 . Near the center of the exhaust. By increases 
in a left-handed (LH) polarization sense, which is the expected signature of an intermediate mode. However a clear 
LH polarized intermediate mode has not developed. The outflow Vz (in (g)) is driven by the J x B force (magnetic 
tension) associated with the decrease in the magnetic field Bz and is linked to the jump in Bz by the Walen relation, 
Yzd — V^n = '^\/ Pu^u/ l^oi^zd/ P — ^zu/ Pu)i whcrc the subscripts "u" and "d" indicate upstream and dowstream 
values of the parameters, respectively 126]. The Walen prediction matches the exhaust velocity very well (dot-dashed 
blue curve in (g)) in the regions where Bz decreases gradually outside of the core of the exhaust. However, contrary 
to the Walen prediction, no further increase in the outflow velocity occurs when \Bz\ decreases sharply to zero in the 
center of the exhaust. This is probably because the magnetic tension in this firehose unstable region is zero. The 
consequence is that the reconnection outflow speed is typically ~ 40% slower than the Alfven speed based on the 
asymptotic magnetic field. 

In order to get a better idea of how slow shocks and rotational waves propagate out from the central exhaust to a 
much larger distance (~ lOOd^) than can be achieved in a reconnection simulation, we compare the results of Fig. [2] 
with those from a quasi-lD Riemann problem in Fig. [3] The angle Obn between the upstream magnetic field and 
the shock normal for this simulation was taken to 75°, which is the approximate value just upstream of the exhaust 
boundary in Fig.|2|f). Data from this simulation was presented earlier [18]. The cuts are very similar to those obtained 
from the cuts across the reconnection exhaust but there are differences. In Fig. [sj^a) the plateau in e at 0.25 is much 
more developed and the region of firehose instability is much smaller than in the reconnection simulation. Unlike the 
results from reconnection, the downstream LH rotational wave in Fig. [3|b) has clearly developed. These rotational 
waves have been identified as the downstream dispersive wavetrains of slow shocks [27 . However, the wavelength of 
the oscillations has been shown to depend on the shock transition thickness, which is not consistent with the dispersive 
wavetrain hypothesis [18 . Finally in the Riemann solution the outflow continues to be linked to the variation in 5^, 
consistent with the Walen relation (see the dot-dashed blue curve in Fig.jsj^g)). Again, this difference is likely because 
e does not fall below the marginal firehose stability threshold in the Riemann simulation. 



III. 1-D SHOCK SOLUTIONS 



To understand the structure of the shocks that appear in the reconnection and Riemann simulations, we evaluate 
the possible transitions using the anisotropic Rankine-Hugoniot jump conditions. These follow from the moment 
integration of the Vlasov equation for a monatomic plasma assuming that the off-diagonal pressure components are 
negligible [H |28l ES] : 
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[...]^ indicates the difference between up- and downstream. The equations are written in the deHoffmann-Teller 
frame where the electric field vanishes, E = Vx'^t — ^t^x = 0. Quantities p, V^, V^, B^^ and Qx are the mass 
density, velocity of the bulk flow in the normal direction (f ), velocity of the bulk flow in the tangential direction {y-z 
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plane), normal component of the magnetic field, tangential components of the magnetic field, and the heat fiux in 
the x-direction. Here P = (P|| + 2P±)/3. Even with a known upstream state, these jump relations still have free 
parameters Sd and Qxd foi" which we will use the measured values from the simulations. The observed Qx does not 
jump across the SS transition of interest so we can discard it. 

The anisotropic MHD equations can be further simplified in the case relevant for magnetic reconnection where the 
normal magnetic field is much smaller than the upstream transverse field B^^^. In this case the normal velocity 
Vx is of the order of Cax = B^/ ^/JlqP and is therefore also small compared with the fast mode wave speed. The fast 
mode can then be eliminated from the equations and to lowest order. Eqs. ([2|-([5| become 
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where the constant mass fiux has been written as F = pVx and has been eliminated from the zero constraint on 
the electric field. 

We first address why a SSS solution does not appear in either the reconnection or Riemann solutions. Since there 
is no reliable analytic model for the downstream value of e, we use the values of e from the simulations to explore the 
possible shock transitions. In Fig. [4]^a), we plot e versus B^ from the right half part of the 75° Riemann simulation 
and a similar simulation with Obn = 30° (black curves). In the 30° data Bz decreases steadily from its upstream 
value to zero. The SS transition in this simulation corresponds to a SSS in which Bzd = [18 . In the 75° data B^ 
decreases down to around 0.67 where e = 0.25 and then decreases further on the plateau where e remains at 0.25. In 
this case, as we have discussed earlier, there is no SSS transition. To understand the reason for these differences, we 
study the SS solutions from Eqs. ([6|-([9|. For = 1 and Bzu = 1 we find the traditional SSS solution in which the 
upstream intermediate Mach number Mju = Vxuy/^opu/Bx is unity and the density compression across the shock 
Rp = p/pu^ downstream transverse fiow Vzd^ and Mach number Mjd are given by 
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where Pu is the upstream value of /3 = 2Pu/ p^B^. Note that the solution is independent of the downstream value of 
£. Mjd depends on Ed only because the intermediate wave speed depends on this parameter. 

The disappearance of the SSS solutions is due to the emergence of a new SS solution from Eqs. ([6|-([9|, an Mju = 
Mjd = 1 transition that prevents the SSS solution from reaching Bzd = [IS]- Eqs. ([6|-([9| yield the jump conditions 
on the density and magnetic field. 
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This 11-SS solution is shown in the red hnes in Fig. [4]^a) as obtained from Eqs. ([2|-([5|. Note that the 11-SS does not 
intersect the measured Bz — e curve for the 30° case at a point other than the upstream state, but does for 75°. If the 
11-SS solution intersects the Bz—s curve the SSS solution is no longer accessible. This is because the presence of the 
11-SS solution at a finite value of Bzd signifies the formation of a minimum in the effective potential il^{By^ Bz) that 
describes the SSS transition and a corresponding maximum at the origin where Bz = By = [19 . The potential in 
this case is similar to that shown in Fig. |5|b). The survival of the SSS solution therefore requires that the minimum 
of Ed as determined from either simulation or observational data fall above the minimum of £ that is accessible by the 
11-SS. This minimum is obtained from Eq. ([l4| by taking Bzd = 0- Thus, the SSS can form as long as the minimum 
of e across the outflow exhaust exceeds 

Thus, in the case of anti-parallel reconnection where counterstreaming ions reduce £ below this critical value there is 
no SSS transition. On the other hand, it is also evident from Fig. |4]^a) that the 11-SS can not produce a transition 



below the minimum value given in Eq. (15), which lies above the 0.25 measured in the simulations. 

The only option remaining for finding a SS solution that can transition to £d = 0.25 is to increase the value of the 
upstream Mach number Mju above unity. That increasing Mju facilitates a SS transition to a lower value of £d and 
is shown in Fig. |4jb), which shows the Bz—e relation obtained from Eqs. (|2|-([5| for Miu = 1.29. The value of M/^ 
was increased until the shock transition curve intersected the black curve at the red circle {Bz^d = 0.67, Ed = 0.25). 
This super-intermediate SS transition, which was earlier denoted as an anomalous slow shock (A-SS) [22] is a possible 
solution for the SS seen in the 75° Riemann simulation. A similar solution with Mju = 1.17 that results in an 
intersection point {Bz^d = 0.65, = 0.25) is shown with the data from the reconnection simulation in Fig. |4]^c). 

To further test the hypothesis that the A-SS describes the reconnection exhaust boundary, we plot the downstream 
predictions of the A-SS transition based on Eqs. ([2|-([5| on top of the data from the Riemann and reconnection 
simulations in Fig. [2] and [3j The predicted jumps of the A-SSs appear as dotted black lines, and they should be 
compared with the solid black lines. Overall, the predicted jumps of these A-SS's agree very well with the observed 
transitions. With the measured inflow speeds Vix^ which are ^ 0.3 and ~ 0.11 in Figs.|2|g) andjsjg), respectively, we 
can calculate the shock speeds of these A-SS's in the rest frame of the far upstream plasma. The local intermediate 
Mach number Mj is then derived and plotted in (h). The choice of Mju being slightly greater than unity and the 
shock therefore being super-intermediate is consistent with both the reconnection and Riemann simulations but is 
most evident in the Riemann case. The measured downstream Mach number in Fig. |2](h), although oscillating, is 
comparable to the prediction, while it agrees very well in Fig. [s] (h) . Our interpretation of why the shock is super- 
intermediate is that the ions stream along the local magnetic field and into the upstream at the Alfven speed while 
encountering a weak infiowing plasma. The interaction between these two oppositely moving constituents would 
naturally drive the formation of a weak super-intermediate slow shock rather than the conventional standing slow 
shock. 

A better understanding of the A-SS transition can be obtained from the pseudo-potential that describes the shock 
transition. This potential describes the variation of as a function of the space variable in the 2-D By — Bz plane. In 
analogy with 2-D particle dynamics, the space variable can be treated as time and the By — Bz plane can be treated as 
a 2-D configuration space. We calculate the potential using the anisotropic derivative nonlinear- Schrodinger-Burgers 
equation [19^. To do this we need an analytic relation between Bz and £. For simplicity we use a straight line from the 
upstream-to downstream states shown in Fig. [sj^b) (Fig. [sja)). The resulting pseudo-potential is shown in Fig. |5|b). 
In this potential a pseudo-particle starts at an upstream state Bzu and slides down the hill to Bzd with By = 0. The 
peak in the potential around Bz = is evident. This potential hill prevents the A-SS from making a transition to 
Bz = 0. The pseudo-particle could, in principle, follow the valley by rotating Bz into By. However, as discussed 
below, this rotation does not take place at this location but further downstream. 

The A-SS shock speed (solid red) and the variation of the intermediate A/ (dot-dashed blue) and slow Xsl (dot- 
dashed red) wave speeds across the transition are shown in Fig. jsjc). All velocities are shown in the upstream 
intermediate frame. Note that upstream the A-SS is super-intermediate and super-slow (the shock speed at Bzu 
exceeds the intermediate and slow characteristic speeds) and that downstream the A-SS is super-intermediate and 
sub-slow (the shock speed at Bzd exceeds the intermediate but is below the slow-characteristic speed). Thus, this 
A-SS is an Mju > 1 to Mjd > 1 transition with Mjd > Mju as is evident in panel (h) of Figs. |2] and [sj 

Two other distinct features of the exhaust boundaries are also evident in Figs. [2]and|3] In the e = 0.25 pedestal 
region the transverse magnetic field Bz continues to decrease and the exhaust velocity continues to increase, the 
latter being consistent with the Walen prediction. Finally, further towards the symmetry line the transverse magnetic 
field Bz abruptly begins to rotate into the y-direction. The pedestal region with constant s = 0.25 can again be 
explored with Eqs. (|6|-([9|. This case is a special case of the 11-SS solution discussed earlier. As in that solution. 
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Ml = ^jV^idQ/epB"^ is a constant so that p (x e is constant. The energy equation then takes the form 
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This equation reveals the unusual conditions that arise when £ = 0.25. This SS solution can support arbitrary changes 
in the transverse magnetic field and associated pressure (determined by the constancy of the magnetic plus plasma 
pressure). The existence of this transition is marked by the short horizontal red segment at 0.25 in Fig. [4]^b). The 
transition is less evident in Fig. |4jc) because of the limited spatial extension in the normal direction of the reconnection 
simulation. Since p is a constant the transverse flow is given by the simplified Walen relation [14] = ^/e/|M)p[BzYd. 
The condition of constant 6 = 0.25 corresponds to the degeneracy of slow and intermediate waves [19]. We therefore 
refer to this solution as a degenerate slow shock (D-SS). Note that an exactly horizontal transition dX Bz — s plane 
{i.e., caused by scattering [18l [19]) can only take place at e = 0.25 in the form of the D-SS. The pseudo-potential 
can again be constructed for the D-SS. The chosen — s relation is shown by the dashed red curve in Fig. (sjd). 
It is a horizontal line from B^u to B^d and then is taken to 6 = at 5^ =0. The latter segment is meant to be 
illustrative. The structure of the pseudo-potential and characteristic speeds shown in Figs. |5|e-f) are insensitive to 
the detailed variation of s as long as it drops below 0.25. Fig. [sje) reveals that the potential is fiat on the e = 0.25 
plateau and again has a barrier around the origin that prevents a direction transition to Bz = 0. On the plateau the 
shock, intermediate and slow wave speeds are equal to the upstream intermediate speed (Fig. [sjf)). The predicted 
weak jumps of D-SS (dotted red) of the Riemann simulation in Fig. [s] are added to jumps of A-SS (dotted black) for 
further comparison with data. The agreement between the simulation data and the D-SS predictions are reasonably 
good. 

The data from both the reconnection and Riemann simulations reveal the remnant upstream tangential magnetic 
Bz abruptly begins a LH rotation at the lower edge of the e plateau. This is especially evident in the Riemann 
simulations. The transition of the coplanar D-SS to this LH rotational mode follows from the structure of the pseudo- 
potential in Fig. p[e). The pseudo-particle moves toward lower Bz but must rotate as a result of the barrier around 
the origin. Fig.]5[f) reveals that, as expected, the downstream intermediate characteristic speed decreases with 
decreasing s and it becomes slower than the slow characteristic speed. The rotational mode is not an intermediate 
shock (IS) which must be super-intermediate upstream and sub-intermediate downstream {i.e., the blue portion in 
Fig. [sjc)). Thus, the rotational wave should be interpreted as a rotational discontinuity (RD) that is linked to the 
D-SS, forming a compound D-SS/RD. The rate of rotation in the Riemann simulation is significantly higher than that 
in the reconnection simulation. This is likely because e falls to zero in the core of the exhaust in the reconnection 
simulation, which reduces the rotation rate of the LH wave. 

The dynamical picture leading to these multiple structures can be understood as following: The counter-streaming 
ions drives the overall tendency to form a A-SS + RD wave (Fig. [sjc)). However, due to rotation-induced scattering 
the e tends to form a plateau at 0.25 [18l |19] and the downstream region locally evolves into a D-SS/RD compound 
wave (Fig. |5|f)). 



IV. IMPLICATION 

For a density of l.Ocm"^, typical of the magnetotail, Idi = 150/cm, which implies that the simulation shown in 
Fig. 1 corresponds to roughly 20Re. During times when the tail has a single large X-line, the signatures discussed 
in this paper should be observable by Geotail and THEMIS. Specifically, we predict that no switch- off- slow shocks 
should develop if the measured values of e across the exhaust fall below the threshold Scr given in Eq. ([T]) . We predict 
that the exhaust boundary will exhibit two distinct SS transitions. The first is an anomalous slow shock (A-SS) that 
is a super-intermediate to super-intermediate transition that drops e from unity to around 0.25. The second is a 
degenerate slow-shock (D-SS) /rotational-discontinuity (RD) compound wave that forms on and maintains a plateau 
with e = 0.25. When e = 0.25 the slow and intermediate wave speeds are equal and the resulting D-SS can support 
an isobaric change in magnetic and particle pressure while leaving the density constant. The RD rotates the residual 
transverse magnetic field in the left-hand sense. Since the cross-exhaust width of the firehose unstable region does not 
increase with distance downstream of the X-line in reconnection simulation, we suggest that the dissipation region 
proper is the source of the firehose unstable region. This firehose unstable region causes the failure of Walen relation. 

Several caveats must be kept in mind. First, because of the relatively narrow angle of the reconnection outfiow 
exhaust, small reconnection events might be too narrow for the s = 0.25 plateau to develop. Second, we have only 
considered 2D simulations that are invariant in the y direction. A more realistic 3D system may develop other features 
that mask or even change our conclusions. Finally, unlike the anti-parallel case, with a significant ambient initial 
guide field By the firehose parameter e does not become as small and as will be shown in a future manuscript a pair 
of RD's switch off the reconnecting magnetic field {Bz). 
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Appendix: Exact solutions of the anisotropic MHD equations 
Equations ([2|-([5| have analytic solutions, 
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where Qx has been neglected. Using the identity 
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with given ^^i, and A^^^, we can plot £ versus Bt and obtain the following jump conditions. 
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In the very oblique limit, we look for solution of A^^ = Su and = e (the 11-SS solution). Bf = is further 
imposed, since this point gives us the minimum Scr of 11-SS transitions. We get 
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Plug these into Eq. (17), choose the plus sign and expand it, the minimum Scr value that a SSS (with a given £„ and 
/?„) can transition to is obtained 



A. 



-2el + (15/?„ 



Eq. is re-derived when £„ = 1, 



|6| 10£„ + 15/3„ + 5 



(30) 



(31) 



We noticed that Eu —Scr > vanishes when Eu approaches value 0.25 from unity. Therefore £ = 0.25 is the solution of 
£u = ^cr- From this we can conclude that a SSS with a > 0.25 and an arbitrary Pu can not transition to £ < 0.25 
region. In other words, 0.25 is an "absolute" barrier across which a SSS with Su > 0.25 can not cross. 
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FIG. 2: From the reconnection simulation profiles of physical quantities along the black cut of Fig.[T]^b). e and Qx in (a), 5's 
in (b), T's in (c), P's in (d), rii and rie in (e), /3 and Obn in (f), Vs in (g) and Mi in (h). Obn is the angle between the 
local magnetic field and the x-direction. The dot-dashed blue curve in (g) is the Walen relation. The dotted black curves are 
the predicted jumps of an anomalous slow shock with Miu = 1.17, and these dotted curves should be compared with the solid 
black curves. 
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FIG. 4: The Bz — s space where the downstream direction is toward Bz = 0. In (a) the lower black curve is measured 
from the right half part of data from the Riemann simulation in Fig. [3] and the upper curve is from a similar simulation with 
Obn — 30°. The red curves show possible transitions of a Miu — Mjd — 1 slow shock transition (11-SS) from the anisotropic 
MHD equations. If the red 11-SS curve intersects the black curve at any location other than in the upstream state, the switch- 
off-slow shock transition is not possible. In (b) the 75° Riemann data in black, the 11-SS transition in red and a Miu — 1.29 
SS transition (A-SS) in blue. The intersection of the black and the blue curves at the red circle {Bz = 0.67, £ = 0.25) is the 
chosen A-SS downstream since this solution allows a transition from e unity upstream to a downstream value of 0.25. In (c) 
plots from the reconnection simulation similar to those in (b) for the Riemann 75° case. The red circle [Bz — 0.65, £ = 0.25) 
is the chosen downstream state. 




FIG. 5: Analysis of the 75° Riemann simulation. In (a) the Bz — £ data from the simulation with a chosen analytic relation 
for £{Bz) linking the upstream and downstream states. In (b) the pseudo-potential ^ from anisotropic derivative nonlinear- 
Schrodinger-Burgers equations [19] showing the potential minimum at the downstream solution and the barrier preventing the 
development of the SSS. In (c) the shock speed (red solid) and the slow (dot-dashed red) and intermediate (dot-dashed blue) 
characteristics, \sl and A/ respectively, show that the A-SS solution is a super-intermediate to super-intermediate transition 
but a super-slow to sub-slow transition. In (d) similar to (a) but with an analytic relation for £{Bz) describing the e = 0.25 
plateau region and extension toward £ = 0. In (e) the pseudo-potential for the D-SS showing that the potential is flat between 
the up and downstream states and the barrier that facilitates a left hand rotation that describes the compound D-SS/RD 
transition. In (f) wave characteristics for the D-SS/RD as shown in (c). 



